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We study spontaneous symmetry breaking in a system of spinless fermions in the Honeycomb 
lattice paying special emphasis to the role of an enlarged unit cell on time reversal symmetry broken 
phases. We use a tight binding model with nearest neighbor hopping t and Hubbard interaction 
Vi and Vb and extract the phase diagram as a function of electron density and interaction within 
a mean field variational approach. The analysis completes the previous work done in Phys. Rev. 
Lett. 107, 106402 (2011) where phases with non-trivial topological properties were found with only 
a nearest neighbor interaction V± in the absence of charge decouplings. We see that the topological 
phases are suppressed by the presence of metallic charge density fluctuations. The addition of next 
to nearest neighbor interaction V2 restores the topological non-trivial phases. 



I. INTRODUCTION 

Topological phases of matter are a new paradigm in condensed matter physics P^l These phases evade a standard 
classification in terms of local order parameters and broken symmetries, being described by topological invariants. 
In addition to the obvious interest from a fundamental viewpoint, the robustness of topological properties against 
certain local perturbations make these phases appealing also in applied physics. The recently discovered time reversal 
invariant topological insulators are a raising star in the family, with topological invariants protected against non 
magnetic disorder by time reversal symmetry (T) P^ Together with time reversal invariant topological insulators, the 
quantum Hall insulating state is the paradigmatic example of a topological phaseP The non trivial topology is in this 
case driven by an external magnetic field which breaks T. In quantum anomalous Hall (QAH) insulators, another 
example of a topological phase, T is broken spontaneously, and a quantized anomalous Hall (AH) conductivity arises 
in the absence of any external magnetic fieldP When the Fermi level does not fall into the bulk band gap there is a 
non quantized contribution to the AH conductivity characterized as a property of the Fermi surface through its Berry 
phase and the systems are termed topological Fermi liquidsP^ 

The Honeycomb lattice is perhaps the best studied case for its special properties! 11 * 12 * It is bipartite and yet due 
to its topology it is easily amenable to frustration. It is fair to say that the pioneer works proposing this lattice as 
a "condensed matter simulation of a three-dimensional anomaly 'EIED together with the analysis of the spin orbit in 
graphene done in refs. H4|15l opened the modern field of topological insulators. In these pioneering works the breaking 
of time reversal symmetry that allowed the topological non-trivial phases was explicitly put in the Hamiltonian in the 
form of complex hopping parameters. The non interacting behavior of the topological phases in the insulating family 
is at present fairl y well understood, and the attention is now shifting to the effect of electron-electron interactions 
on these phasesE3t2H anc [ to the nature of the phase transitions between topological and ordinary phases (see^l and 
references therein). In most cases following the original work by Haldane 2 breaking T is associated with bond order 
(complex hoppings) originating finite flux states! 25 " 29 ! i n this later context the electron spin is not meant to be a key 
ingredient, contrary to the topological insulators known to date where the strong spin-orbit coupling is responsible for 
the non trivial topology. Focusing in the two-dimensional (2D) case, and using spinless fermionic models, the strategy 
then is to search for spontaneously broken T phases showing an AH or QAH effect driven by Coulomb interactions. 

A realization of the Haldane model through electron-electron interactions was obtained at mean field level in the 
honeycomb lattice in refs. I25|29l by adding second neighbor Coulomb interactions. Other proposals involve more 
complex lattices which allow for intracell fluxes, as the checkerboard, the Kagome, or the dec orate d honeycomb 
T-broken superconducting states have also been proposed recently on the honeycomb latticeP^l 

In a recent publication^ 2 W e proposed enlarging the unit cell of simple lattice models as an alternative way to 
drive the spontaneous appearance of phases with broken T. We explored as an example the nearest neighbor (NN) 
tight binding model for spinless fermions interacting through a NN Coulomb interaction in the honeycomb lattice. 
Our main motivation was the idea that enlarging the unit cell - to enable for instance Kekule type of distortions - 
would allow for non-trivial topological phases without the necessity to go to longer range hopping or interactions. 
We also focussed on a region of high doping near the Van Hove singularity where short range electronic interactions 
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are enhanced and give rise to interesting phases. For this purpose we considered a minimal model and restricted the 
mean field decoupling to order parameters of the type < a^bj > ignoring possible charge ordered phases with order 
parameters of the type < cfci >. The result was that T-broken phases with interesting topological features appeared 
at large fillings above the Van Hove filling for reasonable values of the interaction. The most interesting region of the 
phase diagram occurred around the commensurate value n = 1 corresponding to have one electrons per enlarged (six 
atoms) unit cell. There the system was insulating above a critical value of the NN Hubbard interaction V = V±. In 
this work we complete the former analysis by allowing charge decouplings in the mean field equations. We find that 
charge inhomogeneous phases dominate the phase diagram in the region where T broken phases occurred. T broken 
phases reemerge and are stabilized by the addition of a next nearest neighbor (NNN) interaction V^. The results then 
are similar to these obtained at half filling in the pioneer work of of Ref. [25] in the sense that V2 ^ is needed to 
stabilize non-trivial phases. In the present work, however, we have found that the high doping T-broken phases are 
stabilized for V2 <C V\ while at half filling V2 > V\ is required. 

The article is organized as follows: In section [IT] we describe the model and the method of calculation. In Sec. 
[TTT] we explain the phase diagram where for completeness, and to compare with previous results, we present also the 



situation at half filling (III A). In Sec. IIIB we analyze the modification introduced in the V2 = case by the charge 
decoupling at higher fillings. We will see that the charge modulated phases wash out the topologically non-trivial 
phases. Finally we see how these are restored by the inclusion of the second neighbor interaction and describe the 
T broken phases and their band structure. In Sec. |IV| we summarize the situation and discuss some open problems. 
The Appendix [A] contains the technical details of the calculation. 



II. MODEL 



We consider the model describing spinless electrons on the honeycomb lattice, 

where t is the NN hopping, and V\ and V2 the NN and NNN repulsion, all of them real. The operator c\ (q) creates 
(annihilates) a fermion at site z, the number operator is ni = c\ci, and the sums run over either NN sites or 
NNN sites 

We use a 6— atom unit cell to allow for finite flux also in NN loops, and not only in NNN loops. The basis vectors can 
be chosen as ai = 4f (— 1) and a 2 = ^(y/S, 1), and their counterparts in reciprocal space are bi = ^§^(~ 1? V%) 

and b 2 = ^§^(1? V^). The number of Fourier components of local operators is six, aj k = ^= ^2 ieAi ^e 1 ^' ' Vi and 
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b] ^ = ^2 ie B L 4 e2k " Vi j w ^h 6 = 1?2,3. In terms of the new operators the Hamiltonian in Eq. ([I]) reads as 



ff = -*E fl U7k*)\k + H.c. 
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where summation over repeated indices is assumed (thus the factor 1/2 in the last two terms), with z, j = 1, 2,3, and 
where we have defined the 3x3 matrices 
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Note that while a q and /3 q are hermitian, the matrix 7 q is not. 
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III. PHASE DIAGRAM 



The phase diagram of the model is obtained within the variational mean field approach. The details of the mean 
field decoupling and related equations are extensively explained in Appendix A. In brief, we replace the four fermion 
interaction terms in Eq. (2) with bilinears which, when written in real space, can be interpreted as the most general 
hopping and potential energy terms compatible with the reduced translational symmetry of the lattice (6 atom unit 
cell). In total there are 9 NN (£) and 18 NNN (9 per sublattice, \ A an d X B ) complex hopping parameters, which 
are depicted in Fig. 1. In addition there are 6 local energy terms (3 per sublattice, p A and p B ), of which only 5 are 
independent due to charge conservation. See Appendix Al for details. 

Using the variational mean field approach one finds the set of 33 = 3x9+6 mean field equations which, complemented 
by charge conservation, determine the mean field parameters and the chemical potential \i (see Appendix A2 and A3). 
The mean field equations read as, 
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where A^' t *, A^* are phase factors analogous to 7^ defined in Eq (3), r, 
with c = A, B. Detailed expressions for these matrices can be found in the appendix A2l The notation {■■■) MF 
means average in the macrocanonical ensamble taking the mean field Hamiltonian in the Boltzmann factor. 
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Figure 2: (Left) Mean field phase diagram for the half filling case reproduced from ref. [29] The various phases are described 
in the text. SM means semi- metal. (Right) Mean field phase diagram obtained in present work. Lines are guides to the eyes. 
CMs stands for the charge modulated phase discussed in the text. 



In order to obtain the mean field phase diagram we solve the mean field equations self consistently, (see Appendix A3) 
and take the solution (if more than one is obtained) which minimizes the free energy in Eq. ( A1Q| ) (see Appendix A4). 
Care must be taken with charge like order parameters, Eqs. (9) and (10). Due to the frustration introduced by NNN 
interaction, these order parameters may flow to a non-selfconsistent solution where the charge like order parameters 
in different sublattices interchange at each step. Apart from this subtlety, getting a solution is straightforward. 

We will analyze first the phase diagram obtained at half filling which is interesting on its own and later discuss the 
modification introduced in the V2 = case by the charge decoupling. We will see that the charge modulated phases 
wash out the topologically non-trivial phases. Finally we see how these are restored by the inclusion of the second 
neighbor interaction. 



A. Half-filling 



Let us first analyze the half- filled case, where n = + n# — 3 = 0. This case provides a test to the present mean 
field analysis, since a similar approach, also using a 6— atom unit cell, has been taken in ref. [23 For comparison, we 
show the phase diagram obtained in ref. |29]in the left panel of Fig. [2] In the right panel of Fig. [2] we can see the phase 
diagram of the present work (we use the same color code). We plot the different phases (that will be described in 
what follows) as a function of the interaction strength V\ and V2 in units of the hopping parameter t. The half filled 
case was first explored in the original lattice in ref. [25] and non-trivial topological phases were already encountered 
for values of the interaction V2 > V\. 

For V\ < 1M and V2 < 2t the two phase diagrams coincide. For V\ > 1.5t, however, we find that the semi- 
metallic (SM) and the charge density wave (CDW) phases are robust against the Kekule phase. The Kekule phase is 
characterized by an alternating bond strength as shown schematically in the inset of the left hand side of Fig. [3] This 
distortion is important in the physics of graphene because it opens a gap in the spectrum breaking the translational 
symmetry of the original Honeycomb lattice while preserving time reversal (T) and inversion (X) symmetries^. It 
also plays a key role in the models of charge fractionalization in the Honeycomb lattice^. In our approach the Kekule 
phase only appears at much higher V\ and V2. For V2 > 2t a new phase sets in, not predicted in ref. 29. This is a 
charge density wave with reduced rotational symmetry; to distinguish from CDW we denote it as charge modulated 
with modulation also over the sublattice (CMs). In the CDW there is a charge imbalance between sublattices, 
but no inhomogeneity over the sublattice: the charge-like order parameters (pf, pf , p^, pf , p^, pf ) take the form 
(p, —p, p, —p, p, —p). For CMs, however, the charge is modulated also over the sublattice and the charge like order 
parameters take the form (p^—p^p^—p^—p — A,p + A). We note that for V2 > V\ such modulation is naturally 
expected from a classical point of view (remember that the Hartree contribution has a classical interpretation): The 
staggered charge modulation of CDW minimizes the energy coming from NN repulsion, but it does not affect the 
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Figure 3: (Left) Phase diagram with V2 = when charge inhomogeneous phases are not allowed. Legend: (S) symmetric 
phase, i.e. bare graphene with a uniform renormalization of the hopping; (K) Kekule distortion with hopping renormalization 
as shown in the inset; (P) Pomeranchuk distortion of the Fermi surface and hopping renormalization as shown in the inset; 
(K+P) coexistence of Kekule and Pomeranchuk distortions; (T-I) and (T-II) T broken phases discussed at length in the text; 
(RS) broken symmetry state with real hopping parameters, the distortion is neither Kekule type nor Pomeranchuk (reduced 
symmetry). (Right) Phase diagram with V2 = and charge inhomogeneous phases allowed. CM stands for the new charge 
modulated phase discussed in the text. In the CM+RS phase there is also a real, asymmetric renormalization of the hoppings. 

NNN contribution; for V2 > V\ it becomes energetically favorable to reverse a NN dimmer, paying the corresponding 
NN energy ex Vi, but reducing the NNN energetic contribution ex V2. Such inhomogeneous charge modulation over 
the sublattice was not allowed in ref. [29j and that is the reason why the QAH phase dominates over a larger region 
of the phase diagram. 

B. Higher filling 

We will now discuss the results obtained for higher dopings. The doped system has been the subject of atten- 
tion r ecently due to the experimental availability^ and to the interesting phases that emerge near the Van Hove 
fillinjp^IMH] In ref. [35] we have shown that for the Hamiltonian in Eq. ([I]) with V2 = 0, and ignoring charge inho- 
mogeneous phases, T— broken phases show up in the phase diagram V vs n. The phase diagram is shown in Fig. [3] 
(left), where the meaning of the various phases are described. The dominant phases in the region of interest above 
the Van Hove filling (n=0.75 in our units) were a Pomeranchuk instability (P) characterized by rotational symmetry 
breaking as indicated in the inset, and the T— broken phases denoted by TI and Til obtained at 1 < n < 1.5 and 
V\ >2t. These phases are the same as those obtained in ref. 1321 

When charge inhomogeneous phases are allowed we obtain the phase diagram shown in Fig. [3] (right). As can be 
seen, the T— broken phases are washed out by a charge modulated (CM) phase. This charge modulation corresponds 
to a charge imbalance between the two sublattices, and is homogeneous over the sublattice; at half-filling this is 
the CDW discussed in the previous section. This modulation induces a trivial gap between bands 3 and 4 at the T 
point. The system is thus a trivial insulator at n = 0, where the Fermi level falls into the gap, and is metallic for 
n > 0. We see that the CM phase dominates the region of the phase diagram where Pomeranchuk and T— broken 
phases were stabilized before (left hand side of Fig. [3|. The charge decoupling leaves only a small region around the 
Van Hove filling (n=0.75) and for very small values of the interaction where the Pomeranchuk instability is still the 
most favorable phase. At larger values of the interaction (CM+P) the charge modulated phase is accompanied of a 
renormalization of the hopping with the Pomeranchuk symmetry. This is similar to the (CM+RS) phase that occurs 
in the phase diagram around n = 0.75 for larger values of the doping and interaction. 

A natural way to restore the topologically non-trivial phases is to add a NNN interaction V2 • Since we are looking 
for topological phases we have centered our attention in the doping region around the commensurate value n = 1. 

In Fig. [4] we show the phase diagram for V2 > obtained at fixed n = 1.2. Increasing V2 frustrates the CM phase, 
and T— broken phases are recovered. We note that this happens already for V2 <C V\. If V2 is further increased the 
system falls into the CMs phase (charge modulated with modulation also over the sublattice) already encountered at 
half-filling that was discussed in the previous section. 

The new T— broken phases T-I and T-II are similar to the ones described in ref. [32j Fig. [5] shows the real space 
hopping pattern for T-I and T-II respectively. The direction of the arrows represents the sign of the phase of the given 
complex hopping. As discussed in ref. [10] the possibility of having non-trivial topological phases characterized by a 
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Figure 4: Phase diagram in the plane Vi vs Vi at n = 1.2. 



(a) (b) 




Figure 5: (Left) Flux pattern in the T-I phase, (a): The case with V2 = and no charge decoupling in the phase digram at 
the left hand side of Fig. [3] (b) and (c): Flux patterns for the T-I phase in Fig. [4] (Right) Same for the T-II phase. 



finite Hall conductivity in T— broken systems is determined by the discrete symmetries T— and space inversion X— 
preserved in the system. As happened in the case dicussed in ref. |32j the T-II phase is still invariant under inversion 
hence breaking the product T.X and is in principle topologically non-trivial with finite Hall conductivity!^. 

Fig. [6] shows a typical band structure for the T-I and T-II phases. The T-I case is shown in Fig. [6^1eft), where as 
an example we took the point V\ = bt and V2 = 0.7bt. In the right panel of Fig. [6] we show the T-II case, for V\ — 6t 
and V2 = 0.75t. The two figures are very similar to their counterparts when V2 = Cp^. The interesting feature is that, 
unlike what happened there, we do not have any insulating phase away of half filling, even for the commensurate 
value n = 1 whose band structure is very similar to these in Fig. [6j 



IV. CONCLUSIONS AND DISCUSSION 



One of the important points of this and the related work in ref. is the proposal that non trivial topological 
phases can be spontaneously generated from interactions in lattice fermionic systems when translational symmetry is 
relaxed, i.e. when phases with enlarged unit cells are allowed. The employment of enlarged unit cell allows for novel 
intra-cell current patterns emerging from interactions, and thus novel topological phases. We used the Honeycomb 
lattice with extended Hubbard interaction to exemplify the type of new physics expected. In|32j we checked 39 that no 
topological phase can arise in the original lattice (two atoms unit cell), without including a second nearest neighbor 
interactions V2 to allow for the formation of the necessary current loops that break time reversal symmetry. 

We have completed the analysis of ref. |32] including a charge decoupling order parameter. The phase diagram 
is dominated by charge modulated phases with no particular symmetry and no topological phases were found. The 
charge modulated phases in the case V2 = are usual CDW phases but we keep the name CM to emphasize the fact 
that, away from half filling, these phases are metallic. When V2 7^ a different charge modulated phase arise that we 
denote CMs in which the charge is modulated also over the sublattice and the charge-like order parameters take the 
form (p, — p, p, — p, — p — A, p + A). An interesting analysis of charge modulated phases in graphene has been done 
recently irl^. 

In the previous analysis 3 -^ there were two special filling values apart from half filling: The Van Hove filling and 
the commensurate value n = 1 corresponding to four electron per enlarged unit cell. Around the VH density the 
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Figure 6: Mean field band structure in the T-I phase for Vi = 5t and V2 = 0.75£ (left) and in the T-II phase for V\ — 6t and 
V2 = 0.75£ (right). The middle black hexagon indicates the position of the Fermi level. 

system supported a Pomeranchuk deformation and the topological phases were established around the commensurate 
n = 1 filling. The system was gapped only at this particular value. In that case we observed a phase separation in 
the region in the phase space of fillings in between these two values with the two extremes being the stable phases^. 
We have seeing that in the charge modulated phase dominating the phase diagram the n = 1 case is not special any 
more: Even though bands 3 and 4 might be non-trivial now (because of V2), the system at n = 1 is not an insulator 
no "full" gap develops as the Fermi level always crosses bands 4 and 5 (see Fig. [6|. This behavior can be understood 
from the strong coupling point of view due to the frustration induced by the competition between V\ and V2 for the 
ground state. The situation is similar to the quarter filling case in one^ and two dimensions^ where the long-range 
Coulomb interaction enhances the metallic behavior. 

Since this is only a mean field analysis it will be interesting to check the stability of this phases under quantum 
fluctuations. 

Inclusion of an extra coupling to next to nearest neighbors V2, restores the topological phases at fillings around 
n = 1. This model was already studied in the literature, in particular in ref.^ where topological phases were found at 
half filling for values of the interactions V2 > V\ • Allowing the charge decoupling that gives rise to the inhomogeneous 
charge modulation over the sublattice described along this work, favors the T— broken phases that now set at lower 
values of V2. 

The possibility of superconducting phases have been left out in this analysis. Considering superconducting order 
parameters embedded in enlarged unit cells might make novel topological superconducting phases emerge 7 opening 
new routes to realize these, sometimes elusive, phases of matter. Very recent works suggest the formation of an 
eight-atom unit cell spin density wave close to the van Hove filling ( n ~ 0.75 in our notation)^, or topological 
superconducting order exactly at n=0.75 31 . Whether or not these phases dominate close to n ^ 1 or rather a spin 
Hall effect is energetically more stable is still an open question. The T— broken phases described in this work can be 
difficult to observe in graphene since they occur at high values of the interaction and filling but they can probably be 
tested in cold atom experiments with optical latticepMlZ[ 
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Appendix A: Mean field analysis 
1. Mean field decoupling 

The mean field Hamiltonian we propose is 



H M f = Ho 



k 



/ 


Pi 


£ll 


Xf 2 (k) 


e ^12 


(x^i(k))* 




\ 






Pf 




Xf 2 (k) 


e" ikai Si 


(X3 B i(k))* 






(xf 2 (k))* 


6l 




^22 


X^k) 


e -*k-(ai+a 2 )£ 23 






e- ik - a2 Ci* 2 


(Xf 2 (k))* 






^32 


Xf 3 (k) 






x^i(k) 


e ikai 6i 


(X^s(k))* 


^32 




<?33 




V 


£l*3 


X 3 B i(k) 


e ik.(a 1+ a 2 )^* 3 


(X 2 s 3 (k))* 


^33 


pi 


J 



with 



and 



Xf 2 (k) 
X 2 4 3 (k) 
x^i(k) 

Xf 2 (k) 
Xf 3 (k) 
X 3 B i(k) 



X12 



ik-(ai+a 2 ) A,d , ika 2 A,h 
' c A12 e A12 



e -ik-( ai +a 2 )^A,^ 



X23 



A,h 
-X31 > 



B,h 
Xl2 



e ika2 xf 2 rf + e- lk ' ai xf 2 l 



(Al) 



(A2) 



where i^o is the bare Hamiltonian, and we use the spinor notation 1$ 



„ oJ(k) ) 6 t 1 (k),4(k),^(k),4(k),6S(k)]. In 
Fig. [I] we show how the 27 mean field parameters making up the Fock contribution can be interpreted as NN or 
NNN hoppings. The 6 charge like mean field parameters come from the Hartree contribution due to charge imbalance 
between the six sites of the unit cell. Due to charge conservation only 5 of them are independent. 



2. Mean field equations 



Minimizing the free energy functional 

F[£>, Xi p] = (H)mf ~ TSmf = Fmf + (H — H M f) M f — ^mf + (H — H M f) M f + P^e 
we obtain the following set of mean field equations: 



(A3) 



Xij',k 



N 
V 2 



q 



MF 



q \ a j,q a ^q 



MF 



V2 

' N 



q \ h U h ^ 



MF 



V x n B + 3V 2 n A - 3V 2 nf 
Vm A + 3V 2 n B - 3V 2 nf , 



where we have defined ni = ^ J] 



^q \ C i,q C ^q 



and 



It is easy to show, using Eqs. and (A1)-(A2), that the mean field equation for the effective NNN hoppings 

can be cast in the k— independent form 



X % 



V2 
' N 



VA r ' 5 lc ] c- 



MF 
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with 



and 



\A,u _ \A,d _ -iq.(ai+a 2 ) \A,h _ -iq a 2 
A 12,q — 1 A 12,q — C A 12,q — C 



\ A , U _ ^iq-(ai+a 2 ) \ A,d _ i \A,h _ iq-ai 
A 23,q — e A 23,q — 1 A 23,q — e 



A 31,q — e A 31,q — e A 31,q — 1 



\ B >h _ \B,d _ -iq. a2 a#,u _ iq-ai 
A 12,q — 1 A 12,q ~ e A 12,q — e 



\ B >h — J<IZL2 \ B 4 _ a£,u _ iq.(ai+a 2 ) 
A 23,q — e A 23,q — 1 A 23,q — e 



_ p-iq-ai _ i q .(ai+a 2 ) a#,u _ 

A 31,q — e A 31,q ~ e A 31,q ~ 1 



3. Solving the mean field equation 

In order to solve the set of mean field equations given above we need to compute the averages of the form 
&t(q)a^(q)\ . To calculate these averages we introduce the unitary transformation U which diagonalizes HmfQs), 



mf 



UHmf(M)U* = diag[ei(k), . . . , e 6 (k)]. 



The new operators c a (k) are such that 



rci(k)i 




' ai(k) " 


c 2 (k) 




6x(k) 




= 14 




c 5 (k) 




as (k) 


L c 6 (k) J 




L 63 (k) J 



and the average /ftt(q)a i (q)\ and /ct(q)Q(q)\ may be written as 



&J(qK(q) 



MF 



= ($^^a,2j4(q)^C/"| >2i _ 1 c /3 (q) 

\ « /3 
6 

6 

= ^2jt/*,2i-i/[^(q)], 



MF 



MF 



and 



4(q)ai(q) 



6j(q)fti(q) 



a=l 



MF 



MF 



X)^,2i-l^,2i-l/[£a(q)] 
6 
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and the densities (cj(q)ci(q) ) as 

\ / MF 

alfaKfa)) = y] l^,2i-i| 2 /[^(q)] 

\ / MF z — ' 

Q!=l 

[bUdMcO) = £|t/ a , 2 ,| 2 /Mq) 

\ / MF L — ' 



with 

/Mq)] 



exp/3[e a (q) - /i] + 1 
Then we can write the set of mean field equations as, 

6 



& = -^E^E^^a,2i-l/[ea(q)] (A4) 
q a=l 

l ' E A ^ E ^,2i-l^a,2i-l/[£a(q)] (A5) 



A, 5 

Xij 



N 

q a—1 



X*f = -lE^E^^^/^lq)] (A6) 



6 

^ N 

q a=I 

6 

N 



Pt 



V\n B + W 2 n A - ^ E E P a ,2i-i\ 2 /[e a (q)] (A7) 



q a=l 

6 



pf = + 3F 2 n B - ^ E E l^il 2 fle a (d)}, (A8) 



TV 

q a=l 



with 



3 6 



^EEEl^-i| 2 /[ £ «(q)] 

q i=l a—1 

nB = ^EEEl^| 2 /[£ a (q)]. (A9) 

q %—\ a—1 

This set of equations has to be solved self-consistently with the constrain imposed by the Luttinger theorem, which 
reads (ignoring logarithmic corrections in fermion number N e ), 

N e = 1 dn ^ 1 dtt MF = (Af)MF 

N N d/i~ N d/i N 



N 

q,Q! 



and from which we get \i self-consistently. 



4. Free energy 



To check the stability of possible phases we need the Free energy defined in Eq. (A3). For a given set of converged 
order parameters we have, 

T = Q MF + »N e + (H - U M f)mf, (A10) 
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where the mean field grand canonical potential is given by 
and 



q,Q! 



- Umf)mf = V^uaub l^'l 2 - N ^(pt n t + Pf n ?) ~ ( + cx - 



where 



3 



+ |y 2 7Vn A n A - ^A^K* ) 2 + |y 2 JVn B n B - |v 2 7V^(nf ) 2 



TV 



r={A,£} 5={<M,/i} 
r={A,S} 6={u,d,h} 



D 



23 



+ 



D: 



31 



r,5 

12 



+ X23 



23 



r,5 

X31 



31 5 + C ' C -) 



^ = E^q)*^^i-i^, 2i /[^(q)] 

q,ai 

q a = l 

6 

E A tqE^^;,/[^(q)]- 



B,5 



